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ABSTRACT 

We report the results of an X-ray proper motion measurement for the NW 
rim of SN 1006, carried out by comparing Chandra observations from 2001 and 
2012. The NW limb has predominantly thermal X-ray emission, and it is the only 
location in SN 1006 with significant optical emission: a thin, Balmer-dominated 
filament. For most of the NW rim, the proper motion is « 0.30" yr _1 , essentially 
the same as has been measured from the Hcc filament. Isolated regions of the NW 
limb are dominated by nonthermal emission, and here the proper motion is much 
higher, 0.49" yr _1 , close to the value measured in X-rays along the much brighter 
NE limb, where the X-rays are overwhelmingly nonthermal. At the 2.2 kpc 
distance to SN 1006, the proper motions imply shock velocities of ~ 3000 kms” 1 
and ~ 5000 kms -1 in the thermal and nonthermal regions, respectively. A lower 
velocity behind the Ha filament is consistent with the picture that SN 1006 is 
encountering denser gas in the NW, as is also suggested by its overall morphology. 

In the thermally-dominated portion of the X-ray shell, we also see an offset in the 
radial profiles at different energies; the 0.5-0.6keV peak dominated by O VII is 
closer to the shock front than that of the 0.8-3 keV emission — due to the longer 
times for heavier elements to reach ionization states where they produce strong 
X-ray emission. 
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1. Introduction 

The supernova of SN 1006 was observed widely in Asia and Europe, and appears to 
have been the brightest naked-eye supernova ever recorded (Stephenson 2010). Its remnant 
appears in catalog listings as SNR G327.6+14.6, though in this paper we adopt the common 
parlance of using “SN 1006” to refer to the remnant. SN 1006 stands out among historical 
remnants as largest (30'= 19 pc diameter), one of the closest, at 2.2 kpc (Winkler et al. 
2003), farthest above the Galactic plane (550 pc), and least obscured ( N H m 5 x 10 20 cm -2 , 
all of which make it amenable to study in X-rays and other bands. However, it is also the 
faintest of the historical remnants, in total flux and especially in surface brightness, so deep 
observations are required for detailed study. Since SN 1006 is both a prototypical Type la 
SNR, and since it shows some of the best and cleanest evidence for electron acceleration to 
TeV energies in shocks, its close examination can further our understanding of both SNR 
dynamics and particle acceleration in general. 

The overall structure of SN 1006 is that of a nearly circular limb-brightened shell with 
strong bilateral symmetry; its NE and SW limbs are far brighter, in both radio and X-rays, 
than other portions of the rim (Cassam-Chenai et al. 2008; Dyer et al. 2009; Miceli et al. 
2009, and references therein). Koyarna et al. (1995) found that the X-rays from the bright 
NE and SW limbs are not only harder than those from the rest of SN 1006, but also that 
their spectrum is a featureless power-law, a result that provided the first clear evidence 
for diffusive acceleration of charged particles and cemented the long-suspected link between 
supernova remnant (SNR) shocks and cosmic rays (Reynolds 1996). Furthermore, the X-ray 
morphology along the NE and SW rims matches the radio images in exquisite detail (Winkler 
& Long 1997), further confirming a common synchrotron origin for emission in both widely 
separated bands. 

The NW and SE limbs of SN 1006 are dramatically different; the shell in both the 
NW and SE is much fainter and less distinct in the X-ray and radio bands, and the X-ray 


1 Based on observations made with NASA’s Chandra X-ray Observatory, NASA’s Chandra Observatory is 
operated by Smithsonian Astrophysical Observatory under contract # NAS83060 and the data were obtained 
through program GOl-12115. 
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spectrum is soft and thermal, dominated by lines from He-like and H-likc oxygen between 
0.5 and 0.8 keV (Long et al. 2003). The NW limb is unique in that only here can significant 
optical emission be seen: a set of delicate filaments along the outer NW rim, with emission 
that consists solely of hydrogen Balrner lines (van den Bergh 1976; Schweizer & Lasker 1978; 
Ghavamian et al. 2002). 2 First observed in the remnant of Tycho’s supernova (Kirshner 
& Chevalier 1978), B aimer-dominated, or “nonradiative” filaments arise when a fast shock 
expands into a low-density, partially neutral environment. The Balrner emission has two 
components: neutral ff atoms that enter the shock can be excited, and subsequently decay 
to produce narrow lines (with a width characteristic of the pre-shock temperature); or they 
can undergo charge exchange with hot protons to produce broad lines whose width is closely 
related to the post-shock proton temperature (e.g., Chevalier et al. 1980; Ghavamian et al. 
2002; Heng 2010). Since the lifetime of neutral atoms in the hot post-shock environment 
is very short, the Balrner filaments can occur only immediately behind the shock, and thus 
delineate the current position of the shock. 

Proper motions of the brightest Balrner filaments in the NW were first measured by 
ffesser & van den Bergh (1981) using photographic plates taken five years apart, and more 
recently by Long et al. (1988) and by Winkler et al. (2003). The latter measurement, using 
CCD images taken over an 11-year baseline, obtained Una = 0.280 ± 0.008" yr _1 . Based on 
the Balrner profiles and models for the shock conditions, Ghavamian et al. (2002) obtained a 
value of v s = 2890 ± 100 kins'” 1 for the shock velocity in the NW, and combining this value 
with the proper-motion measurement yields a geometrically determined distance of 2.2 ±0.1 
kpc to SN 1006. Having a relatively precise distance of course means that proper motions 
measured by any means, anywhere around the SN 1006 shell, translate directly into a shock 
velocity at that point, which is significant because v s is the most crucial parameter in the 
theory of shock acceleration (e.g., Reynolds 2008). 

A useful way of expressing the proper motion is through the expansion parameter m: 
the power-law index in R oc t m , where t is the age of the remnant. This parameter can be 
interpreted as the ratio of current expansion rate divided by the mean rate over the remnant’s 
lifetime, m = /it/ 9, where 6 is the angular radius. The expansion index gives some insight 
into a remnant’s evolutionary state, e.g., m — 1 for free expansion, or m = 0.4 for Sedov 
expansion. 

Proper-motion measurements of SN 1006 have been made in various bands with different 
precision. Radio measurements of the global expansion gave m ra di 0 = 0.48 ± 0.13 (Moffett 


2 In addition to the NW filaments, there is exceedingly faint and more diffuse Balrner emission surrounding 
almost all of the SN 1006 shell (Winkler & Long 1997; Winkler et al. 2003). 
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et al. 1993), considerably larger than the value of m^a = 0.34 ± 0.01 for the NW optical 
filament alone (Winkler et al. 2003). Since the stronger optical emission in the NW indicates 
a higher external density, it is not surprising that the expansion should be slower than the 
average around the shell, since the higher the density, the more rapidly the expansion will 
decelerate. (If one assumes pressure equilibrium, which may pertain in SN 1006 but is not 
assured, the preshock density will scale inversely as the square root of the shock velocity.) 
In a recent paper, we obtained an X-ray measurement of the proper motions along the E-NE 
limb of SN 1006, using a pair of Chandra ACIS-S observations from 2000 and 2008, and found 
H = 0.48 ± 0.04" yr -1 , or m = 0.54 ±0.05, and v s = 5000 ±400 kms' 1 (Katsuda et al. 2009). 

In this paper, we present the first X-ray measurement of the proper motion along the 
NW limb of SN 1006, using a pair of deep Chandra ACIS-S observations from 2001 and 
2012. These new measurements cover much the same region as the optical proper motion 
measurement by Winkler et al. (2003), and the results are consistent. In addition to the 
thermal X-ray emission that dominates most of the NW limb, there are two small regions 
with a much harder X-ray spectrum, free of emission lines and consistent with nonthermal 
synchrotron emission. We have also measured the motions in these regions and find these 
to be consistent with our earlier measurements (Katsuda et al. 2009) for the much brighter 
synchrotron-dominated NE rim. The 2012 April observation of the NW limb represents the 
first in a series of deep ACIS images covering the entirety SN 1006, being carried out as a 
Chandra Large Project. Further results will be reported in subsequent publications. 


2. Observations and Results 

The initial Chandra observation of the NW limb was carried out with the ACIS-S in 
2001 April and reported by Long et al. (2003). In order to measure proper motions in the 
NW, we carried out a repeat observation in 2012 April; the roll angle and exposure time were 
almost identical for both observations, and the positions were similar, with the second-epoch 
one displaced slightly to the NW in order to include more of the pre-shock background. The 
two observations are summarized in Table 1. In both cases we reprocessed the level- 1 event 
hies with CIAO ver. 4.4 and CALDB ver.4.5.1. 

A three-color image of a portion of the NW limb of SN 1006 as seen with Chandra in 
2012 is shown in Figure 1, where red, green, and blue correspond to 0.5-1 keV, 1-2 keV, and 
2-7 keV bands, respectively. Most of the X-ray emission in the NW is thermal, and appears 
reddish in this image, but there are also a few fainter nonthermal, synchrotron- dominated 
regions that appear more bluish. The difference in spectral character between the thermal 
and nonthermal regions is most apparent in spectra from these regions, shown in Fig. 2. 
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For the sharp, bright shock front in the NE, an outward displacement was evident in 
a simple difference image between the images at epochs eight years apart (Katsuda et al. 
2009). In the NW, however, the surface brightness is much lower and the shock front less 
sharp, so the expansion is less obvious. In order to measure it, we have integrated along the 
shock front to give radial profiles in three thermal regions where the shock is well defined, 
labeled Tl, T2, and T3, and also in two nonthermal regions, labeled NT1 and NT2. 

As described in Sections 2.2 and 2.3, we have measured the proper motions in different 
regions in different energy bands. In order to produce flux-calibrated images, it was necessary 
to consider the variation in effective exposure time with both energy and position on the 
detector. We created exposure maps at monochromatic energies: 0.6 keV, 0.55 keV, 1.2 keV, 
and 2keV for 0.5-1 keV, 0.5-0.6keV, 0.8-3 keV, and 1-8 keV bands, respectively — energies 
that roughly correspond to mean photon energies in the energy bands of interest, given the 
spectral shape in each region (Fig. 2). 


2.1. Image Registration 

Before measuring the proper motion, we first registered the images using as fiducials four 
bright point sources, all with significance levels greater than 10 cr and point-spread-function 
(PSF) sizes less than 2" in both images. 3 We used the CIAO tool wavdetect to obtain 
positions for the sources at the two epochs, which we summarize in Table 2; the sources are 
also indicated in the cyan circles on Fig. 1. Sources PI and P3 correspond to faint ( R to, 19) 
star-like objects in the NOMAD catalog, (Zacharias et al. 2005), both with negligible proper 
motion, while the P2 and P4 sources seem to have no optical companion. 4 All four sources 
have a hard spectrum, and all are likely to be background AGNs; we assume that all have 
negligible proper motions. The offset between the two images, obtained from a simple mean 
of the position differences listed in the two right-hand columns of Table 2, are 0.06" and 
0.18", in R.A. and Deck, respectively. Such offsets are consistent with the stated absolute 
astrometric accuracy for Chandra and the ACIS-S array. 5 After correcting for these offsets, 
the RMS residuals in position difference between the two epochs for the four point sources 
are 0.37" in R.A. and 0.27" in Deck In our proper motion analysis that follows, we have 


3 Here PSF size is defined as a radius in which 39% of total source photons are enclosed. 

4 Source P2 is located 2.5" NW of another faint object in the NOMAD catalog, but since the position dif- 
ference is well outside the error ellipses from both ACIS observations, this is unlikely to be a true association. 
In any case, this object also has no proper motion in the NOMAD catalog. 

5 http:/ / cxc.harvard.edu/cal/ASPECT/celmon/. 
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applied the mean offsets, and have included both statistical uncertainty and uncertainty 

in the registration (for the radial direction), calculated as \J (0.37" sin 6) 2 + (0.27" cos#) 2 , 
where 9 is the azimuth angle for each area, measured counterclockwise from North. 


2.2. Thermal Regions 

We selected regions Tl, T2, and T3 both because they coincide with regions studied 
previously and because the shock front is sharply defined there. Region Tl includes a small, 
bright X-ray knot, located within an Ha bubble ahead of most of the NW shock, and also 
includes the bright X-ray knot where Vink et al. (2003) directly measured both electron and 
ion temperatures using XMM-Newton reflection grating spectra. They demonstrated that T e 
remains much lower than T) well behind the shock, and also that non-equilibrium ionization 
conditions prevail. Regions T2 and T3 trace the shock front delineating the NW limb where 
proper motions of the Ha filaments have been precisely measured (Long et al. 1988; Winkler 
et al. 2003). 

In Figure 3, we plot the projected one-dimensional profiles for epochs 2001 and 2012, 
which clearly show outward motion in all three regions. The profiles have been generated 
from vignetting-corrected images in the 0.5-1 keV band, which contains a large fraction of the 
total X-ray emission. To quantitatively measure the shifts between epochs, we have followed 
the approach we used previously (e.g., Winkler et al. 2003; Katsuda et al. 2009). Briefly, the 
method is to minimize the y 2 value for the difference between the second-epoch profile and 
the shifted first-epoch one, as a function of the amount of the shift. When calculating the y 2 
values, we use data points within the vertical dashed lines indicated in each panel of Figure 3, 
so that we can concentrate on features of interest. The results of these measurements for all 
three regions are given in Table 3. 

The X-ray proper motion we measure in region T2 is consistent, within the uncertainty, 
with the optical measurement for the same region (0.280" ± 0.008" yr -1 : Winkler et al. 2003). 
The proper motions in regions Tl and T3 are slightly higher than in T2, but all are con- 
sistent within the uncertainties. Taking the distance to SN 1006 as 2.2 kpc (Winkler et al. 
2003), the expansion velocity for Tl is 3300 ± 200 ± 300 kms -1 (statistical and registra- 
tion uncertainties, respectively), in reasonable agreement with the velocity inferred from the 
X-ray line width, where Vink et al. (2003) found that the velocity is > 4000 kms -1 , if no 
significant temperature equilibration has taken place, but could be as low as ~ 3000 kms -1 
with some adiabatic cooling. 

The first radio measurement of the SN 1006 expansion was by Moffett et al. (1993), who 
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measured a global average expansion rate of 0.44"±0.13"yr _1 over a baseline of only 8 years. 
They also give local rates in four broad azimuthal sectors, but the NW rim was too indistinct 
to give any measurement at all there. However, Moffett et al. (2004) have since obtained a 
third-epoch radio image of SN 1006, and with a longer baseline obtained expansion rates in 
the NW that are consistent with the optical values. 

In Table 3, we also give the local values for the expansion index, m. Calculating m = 
[it/ 6 requires an angular radius, 6 , and to define this we have taken the expansion center 
to be R.A. (J2000.) = 15:02:54.9, Deck (J2000.) = -41:56:08.9, as determined by Katsuda 
et al. (2009) from the A0514THRI X-ray image mosaic of the entire remnant. 

For Area T2, we have also examined radial profiles and proper motions in distinct energy 
bands. Of the three thermal regions, T2 is the best suited for such a study because it has 
the sharpest shock front and the best photon statistics. As shown in Figure 4 (left), the 
radial profiles vary with X-ray energy; the peak of the 0.5-0. 6 keV (a band dominated by 
O VII triplets) profile is located just behind the Ha: filament, whereas the 0.8-3 keV peak 
is broader and shifted by ~ 15" further behind the shock. This must be primarily due to 
the evolution of the ionization to produce He-like species of elements heavier than oxygen, 
as shown in the profiles and spectra along the NW rim that appear in Long et al. (2003). 
Ionization of higher-Z metals, as well as the gradual increase in the electron temperature 
behind the shock, both make the X-ray spectrum harder toward the center of the SNR. We 
have measured proper motions in these two energy bands, using the one-dimensional profiles 
shown in Figure 4 (center and right), and find no significant energy dependence; the values 
are given in Table 3. 


2.3. Nonthermal Regions 

We have also carried out proper-motion measurements for the small regions in the NW 
that are dominated by nonthermal X-rays, i.e., areas NT1 and NT2 in Figure 1. Since the 
X-ray emission in these areas is harder, we have used a higher and broader energy band, 
1-8 keV, to measure the proper motion. The radial profiles in Figure 5 show clear outward 
motions that are significantly larger than those found for the thermal regions. Measurements 
using the same y 2 minimization technique yield almost identical values of 0.48 to 0.49" yr^ 1 
for both NT1 and NT2, both ~ 50% higher than the values obtained in the thermal regions. 
Furthermore, the proper motions in NW regions NT1 and NT2 are consistent with our 
previous measurements of 0.48 ± 0.04" yr _1 measured at more than a dozen locations along 
the shock in the much brighter synchrotron-dominated NE region of SN 1006 (Katsuda et al. 
2009). 



3. Discussion 


We have confirmed that for the thermal-dominated NW shock, the proper motions 
measured in X-rays are consistent with those measured optically for the Balmer-dominated 
filaments that delineate the shock. While not unexpected, this is the first time that proper 
motions have been measured in both X-ray and optical bands for the same region in SN 1006. 
This is noteworthy because for other young remnants, past measurements in different bands 
have sometimes produced widely discrepant results (e.g., see the summary in Katsuda et al. 
2009), though these may well have been due to comparing local measurements with global 
averages. For the Tycho SNR, at least, recent high- resolution measurements with Chandra 
by Katsuda et al. (2010) are in reasonable agreement with radio measurements by Reynoso 
et al. (1997) around the entire shell. We have now seen that for SN 1006, X-ray and optical 
measurements along the NW rim are in excellent agreement. But as the capability for 
detailed measurements becomes available, it is important to keep in mind that the features 
themselves can evolve, and that new features can appear while others disappear, giving the 
illusion of proper motions in snapshot images taken many years apart. 

Our confirmation that the shock velocity in the thermal-dominated NW limb is slower 
than at other locations around the SNR shell provides direct evidence that the SNR shock 
is interacting with a denser ambient medium in the NW region, as has previously been 
suggested by several studies. Heng et al. (2007) found that the Balmer-dominated Ha 
filaments in the NW (which have significantly higher surface brightness than optical emission 
elsewhere in SN 1006, Winkler et al. 2003) are consistent with a pre-shock density of 0.15-0.3 
cm -3 , higher by a factor of several from that around the rest of the remnant. This density 
is similar to that inferred from X-ray spectra from Chandra (Long et al. 2003) and XMM- 
Newton (Acero et al. 2007; Miceli et al. 2009), as well as from FUSE UV spectra (Korreck 
et al. 2004). Most recently, Winkler et al. (2012) report the detection of 24/jrn infrared 
emission arising from warm dust grains in SN 1006. This emission is seen only along the NW 
rim, and Spitzer spectra of the warm dust is consistent with a post-shock density of 1 cm -3 . 

The shock in the NW of SN 1006 is apparently encountering a localized region of higher 
density than that surrounding the rest of the remnant, possibly representing the start of 
an encounter with the extended H I concentration just beyond the rim found by Dubner 
et al. (2002). As expected where the pre-shock density is higher than average, the rim in 
the NW is indeed flatter than elsewhere, with a radius up to 15% smaller there than the 
average around SN 1006. Since the NW filament lies at a radius only slightly smaller than 
the remnant mean, this slower expansion cannot have been going on for very long, and so 
the shock must have encountered the denser region fairly recently. 

Another point of view is to consider the expansion index m. In the E and NE, m ~ 0.54, 
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indicating that the expansion is well short of reaching the Sedov phase (m = 0.4) in these 
regions, consistent with expansion into a very low-density environment (Katsuda et ah 2009). 
But in the NW, where the radius is only slightly smaller than elsewhere, m ~ 0.38, consistent 
with Sedov expansion, and thus with a higher density. 

The small sections of the NW limb where nonthermal X-ray emission dominates have 
proper motions, and hence shock velocities, significantly higher than those in the thermal 
regions. As noted above, the measured values of ~ 0.5" yr -1 for the nonthermal regions 
NT1 and NT2 give shock velocities ~ 5000 kms -1 at a distance of 2.2 kpc, similar to those 
all along the nonthermal, synchrotron-dominated NE limb (Katsuda et al. 2009). Thus it 
appears likely that in SN 1006, nonthermal-dominated regions are commonly associated with 
fast shocks of ~ 5000 kms -1 . Higher shock speeds than in the thermal regions indicate a 
lower pre-shock density for the the nonthermal regions, suggesting the existence of low- 
density pockets within the generally higher ISM density to the NW, or other small scale 
inhomogeneities. 

Higher shock speeds in the nonthermal regions are also consistent with the theoretical 
view that faster shocks can enhance synchrotron X-ray emission by boosting the roll-off 
frequency zy-oiioff- Diffusive shock acceleration theory predicts a power-law electron energy 
distribution with an exponential cutoff at some characteristic energy, and z/ ro u 0 g is the peak 
frequency emitted by electrons with that energy (e.g., Ellison et al. 2000). Reynolds (2008) 
has shown that v m \\oS oc v* when the maximum energy of accelerated particles is limited by 
the SNR age, or v T0 \\ 0 s oc v% when synchrotron losses limit the maximum energy. 

Still, a shock speed of ~ 3000 kms -1 should be high enough for substantial particle ac- 
celeration to X-ray-emitting energies. The strength and orientation of the upstream magnetic 
held, as well as the electron diffusion coefficient, also play important roles in determining 
acceleration efficiency, and could inhibit particle acceleration in the NW while enhancing it 
in the NE and SW. In addition, a partially neutral pre-shock medium, a condition that the 
existence of bright Balmer filaments indicate must pertain along most of the NW limb, may 
inhibit shock acceleration (Raymond et al. 2011; Blasi et al. 2012). In particular, Blasi et al. 
find that for shock speeds below 3000 kms -1 , the energy spectrum for accelerated particles 
can steepen markedly at low energies, reducing their X-ray synchrotron emission compared 
to that from faster shocks. 

Finally, we note that the SE portion of the remnant has predominantly thermal X-ray 
emission (Rothenhug et al. 2004), but resembles the nonthermal regions in terms of the 
observed local remnant radius (which implies the mean shock velocity over the lifetime of 
the remnant) and in the faintness of its Ha emission (which implies the ambient density). 
Cassam-Chenai et al. (2008) have attributed the dramatic difference largely to a far lower 
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efficiency for particle acceleration to the SE than in the NE or SW. Measuring proper 
motions along the SE rim should be an important diagnostic for understanding the nature 
of particle acceleration in SN 1006. 


4. Summary 

We have measured X-ray proper motions in the NW limb of SN 1006, using the first of 
our new Chandra Large Project observations. For the thermally-dominated X-ray regions 
along most of the NW limb, we find proper motions of 0.28-0. 35" yr _1 , consistent (within the 
uncertainties) with optical measurements for the Balmer filaments, 0.28" yr -1 , by Winkler 
et al. (2003). Even along the thermal-dominated NW rim, however, we find two small, 
isolated regions with harder, non-thermal X-ray emission. Proper motion measurements for 
these regions are much higher: nx ~ 0.49" yr -1 , almost identical to the values measured 
by Katsuda et al. (2009) along the synchrotron-dominated E to NE rim of SN 1006. We 
attribute the slower motion along most of the NW rim to the relatively recent encounter of 
the SNR shock with a higher density pre-shock sheet. However, the existence of small regions 
with non-thermal X-ray spectra and much higher shock speeds in close proximity to the more 
extensive regions with thermal spectra and slower shock speeds suggests low-density pockets 
within this denser sheet, or other small-scale inhomogeneities. 
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Tabic 1. Chandra ACIS-S Observations of SN 1006, NW 


ObsID 

RA(J2000.) 

Dec(J2000.) 

Roll Obs. Date 

Exposure 

PI 

1959 

15:02:22.5 

-41:47:08 

30.2° 2001 Apr 26 

89.0 

Long 

13737 

15:02:15.9 

-41:46:10 

31.7° 2012 Apr 20 

87.1 

Winkler 



Table 2. Positions of Fiducial X-ray Point Sources 
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Table 3. 


Summary of Proper-Motion Measurements 


Region 

Energy Band 
(keV) 

Radius 

0 

Proper Motion 3, 

("yr- 1 ) 

Expansion Index 3 

Velocity 3,6 
( kms” 1 ) 

(X^)min C 

T1 

0.5-1 

14.1 

0.32±0.03±0.03 

0.38±0.03±0.04 

3300±200±300 

0.99 

T2 

0.5-1 

13.3 

0.28±0.02±0.03 

0.34±0.03±0.04 

2800±200±300 

0.58 

T2b 

0.5-0. 6 

13.3 

0.28±0.03±0.03 

0.35±0.04±0.04 

2900±300±300 

0.76 

T2c 

0.8-3 

13.3 

0.28±0.06±0.03 

0.35±0.08±0.04 

2900±600±300 

0.94 

T3 

0.5-1 

14 

0.35±0.05±0.03 

0.42±0.06±0.03 

3700±500±300 

0.84 

NT1 

1-8 

14.7 

0.48±0.05±0.03 

0.54±0.06±0.04 

5000±500±300 

1.32 

NT2 

1-8 

15.2 

0.49±0.04±0.03 

0.53±0.05±0.03 

5000±400±200 

1.13 


a The first and second error terms represent 1 -a statistical and registration uncertainties, respectively. 
b Assumes a distance of 2.2 kpc. 
c Minimum y 2 per degree of freedom. 
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Fig. 1. — Chandra three-color image of the northwestern limb of SN 1006. Red, green, 
and blue correspond to fluxes in the 0.5-1 keV, 1-2 keV, and 2-7 keV bands, respectively. 
The image is binned by 1" and has been smoothed by a Gaussian kernel of a = 6". The 
intensity is scaled as the square root of the count rate. The fields of view of the Chandra 
observations of the northwestern limb (ObsIDs 1959 and 13737 covering relatively inner and 
outer regions, respectively) are shown as white boxes with ACIS chip identifications, and the 
four fiducial sources (PI - P4 from west to east) used to register the images are indicated by 
the cyan circles . The regions where we extract radial profiles are shown in white (thermal 
dominated) and yellow (nonthermal dominated); note the spectral difference between the 
prominent features in the thermal vs nonthermal regions. The field measures 17' x 13' and 
is oriented north up, east left. 
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Fig. 2. — Chandra ACIS spectra extracted from the thermal and non-thermal regions indi- 
cated in Fig. 1. For clarity, spectra from different regions have been scaled by the factors 
indicated. The spectra have been extracted not the entire the regions shown in Fig. 1, but 
instead from the filamentary (T2, T3, and NT2) or knotty (T1 and NT1) features that we 
used for the proper-motion measurement. 
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Fig. 3. — Radial profiles extracted from the thermally-dominated areas shown in Figure 1. 
Data points in black and red represent the 2001 and 2012 epochs, respectively. The intensity 
in 2012 is scaled to equalize that in 2001 by factors of 0.93, 1.01, and 1.11 for Tl, T2, and T3, 
respectively. The dashed vertical lines demarcate the regions used in the y 2 measurement of 
the proper motions. The profile labeled Tl (left panel) covers both the Tl and NT1 regions 
in Figure 1; the smaller peak represents emission from the NT1 region, which is expanding 
faster than the thermal-dominated Tl peak. Figure 5 (left) shows the profile from exactly 
this same area, but for higher energy X-rays. 




Fig. 4. — Left: Radial profiles in the T2 region. Black, red, and blue represent Hcc, 0.5- 
0.6 keV, and 0.8-3 keV emission, respectively. The X-ray data taken in the two epochs are 
combined after taking account of proper motions which amounts to 3 arcsec between the 
two observations. The Ha profile is shifted by an expected proper motion (Winkler et ah 
2003) at the X-ray observation. Center and right: Same as Figure 3, but for different energy 
bands. The intensity in 2012 is scaled to equalize that in 2001 by factors of 1.01 and 0.93 
for 0.5-0.6keV and 0.8-3 keV, respectively. 
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Fig. 5. — Same as Figure 3, but for the NT1 and NT2 regions, using the 1-8 keV band. The 
intensity of 2012 is scaled to equalize that of 2001 by factors of 1.09 and 0.94 for NT1 and 
NT2, respectively. The profile labeled NT1 (left panel) is from exactly the same area as the 
one labeled T1 in Figure 1; above 1 keV, emission from the outer, non-thermal region (NT1) 
is stronger than that from the inner, thermal one (Tl). 


